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^ , ABSTRACT 

■ This paper explores the X-ray properties of 'normal' galaxies using a shallow (2- 

lOkscc) XMM-Newton survey covering an area of « 1.5 deg 2 . The X-ray survey over- 
y—{ ■ laps with the 2dF Galaxy Redshift Survey. Compared with previous studies this has 

J> ' the advantage that high quality spectra and spectral classifications (early, late type) 

exist for all galaxies to bj = 19.4. Moreover, sources with optical spectra revealing 
powerful AGNs can easily be discarded from the normal galaxy sample used here. In 
particular, we present stacking analysis results for about 200 galaxies from the 2dF 
JiJ . Galaxy Redshift Survey at a mean redshift of z » 0.1. We detect a strong signal 

1 for the whole sample (« 6cr) in the soft 0.5-2 keV band corresponding to a flux of 

\ w 7 x 10 -16 erg sec -1 cm -2 and a luminosity of « 2 x 10 40 ergsec^ 1 . A statistically 

— . . significant signal is also detected for both the early and late galaxy sub-samples with 

' X-ray luminosities of « 3 x 10 40 and w 5 x 10 39 erg sec -1 respectively. In contrast, 

P-i, no signal is detected in the hard 2-8 keV band for any of the above samples. The 

Q ■ mean Lx /Lb ratio of the spiral galaxy sample is found to be consistent with both 

| local (< 100 Mpc) and distant (z k, 1) samples (after accounting for differences in 

c/3 . Lb) suggesting little or no evolution of the X-ray emission mechanisms relative to the 

optical. The 0.5-2 keV XRB contribution of the spiral galaxy sub-sample at z ss 0.1 is 
\ estimated to be 0.4 per cent in broad agreement with the XRB fractions estimated in 

• *h . previous studies. Assuming that star-forming galaxies evolve with redshift as (1 + z) k 

' the present data combined with previous studies suggest k < 3. The k values are 

constrained by the relatively low fraction of the soft X-ray background that remains 
unresolved by deep surveys (6-26%). Higher k values will result in an overproduction 
of the soft X-ray background given the fraction already attributed to AGN, groups 
and clusters. The mean X-ray emissivity of spiral galaxies at z ss 0.1 is also estimated 
and is found to be consistent within the uncertainties with that of local H II galaxy 
samples. Using the mean X-ray emissivity of the spiral galaxy sub-sample we esti- 
mate a global star-formation density of 0.009 ± 0.007M Q yr -1 at z « 0.1. Although 
the uncertainty is large this is lower than previous results based on galaxy samples 
selected at different wavelengths. Nevertheless, given the large uncertainties involved 
in converting X-ray luminosity to star-formation rate the agreement is surprising. 

Key words: Surveys - Galaxies: normal - X-rays: galaxies - X-ray:general 



1 INTRODUCTION 

The diffuse X-ray background (XRB), first discovered by Gi- 
acconi et al. (1962), is believed to be dominated by obscured 
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and unobscured AGNs and QSOs. Indeed, deep ROSAT sur- 
veys to the limit ~ 10~ 15 erg s _1 cm -2 have resolved 60- 
80 per cent of the soft band (0.5-2 keV) XRB into discrete 
sources (e.g. Boyle et al. 1994; Hasinger et al. 1998; Lehmann 
et al. 2001). Follow-up spectroscopic studies have demon- 
strated that most of them are AGNs with a small frac- 
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tion associated with groups or clusters (Georgantopoulos 
et al. 1996; Schmidt et al. 1998; Lehmann et al. 2001). Re- 
cently, sensitive Chandra surveys extended the studies above 
to fainter flux limits, ~ 10~ 17 — 10~ 16 ergs -1 cm~ 2 , resolv- 
ing 80-95 per cent of the soft band XRB (Mushotzky et al. 
2000; Giacconi et al. 2002; Rosati et al. 2002). A key de- 
velopment from these ultra-deep surveys is the emergence 
of X-ray faint sources with properties similar to those of 
'normal' galaxies (i.e. systems not obviously dominated by 
AGN activity). These sources, although too faint to domi- 
nate the XRB (< 40 per cent; Giacconi & Zamorani 1987; 
Griffiths & Padovani 1990; Pearson et al. 1997; Alexander 
et al. 2002; Hornschemeier et al. 2002a; Persic & Rephaeli 
2003; Ranalli, Comastri & Setti 2003) are likely to outnum- 
ber AGNs at faint flux limits (< 10 -17 erg sec -1 cm~ 2 ; Horn- 
schemeier et al. 2002b; Ranalli et al. 2003). Study of this 
population to high redshifts provides a census of the global 
star-formation rate complementary to that obtained from 
other wavelengths. More importantly, the X-rays are unique 
in that they are the only tool available to probe the evo- 
lution of low- and high-mass X-ray binaries as well as hot 
gas. Comparison of the X-ray properties of local 'normal' 
galaxy samples with high redshift ones provides a picture of 
the evolution of the Universe that is complementary to that 
inferred from the study of QSOs that dominate the XRB. 

Study of 'normal' galaxies in the local Universe (< 
100 Mpc) has been initiated by both the Einstein and the 
ROSAT satellites. Ellipticals were found to retain large 
amounts (w 10 8 -10 n M Q ) of hot gas (T « 10 7 K) that dom- 
inate their X-ray emission (Fabbiano 1989; Kim, Fabbiano & 
Trinchieri 1992; Fabbiano & Shapley 2002), while the spiral 
galaxy X-rays were shown to originate from a combination of 
stellar point sources and diffuse hot (T w 1 — 8 x 10 6 K) gas 
(Read, Ponman & Strickland 1997; Read & Ponman 2001). 
Moreover, a number of studies have demonstrated the sig- 
nificance of the environment, the star- formation history and 
the total galaxy mass in understanding the 'normal' galaxy 
X-ray emission (Read & Ponman 2001; O'Sullivan, Forbes 
& Ponman 2001; Helsdon et al. 2001; Fabbiano & Shapley 
2002). 

In the last few years the launch of XMM-Newton and 
Chandra, with their significantly improved sensitivity and 
resolution, enabled study of the X-ray properties of 'normal' 
galaxies at cosmological distances. Brandt et al. (2001b) 
studied the mean X-ray properties of bright spirals at z « 
0.5 using a 500 ks Chandra exposure of the HDF-North. 
They found an increase in the mean X-ray luminosity of 
their sample compared to local spirals suggesting evolution. 
Hornschemeier et al. (2002a) used stacking analysis to inves- 
tigate the mean properties of z « 1 spirals in the Chandra 
Deep Field North. They estimated a mean 0.5-2 keV flux 
of « 5 x 10~ 18 ergcm~ 2 s _1 for the spirals in their sample 
and found evidence for an increase in their Lx/Lb ratio. 
However, statistical uncertainties did not allow firm conclu- 
sions to be drawn. Stacking of Lyman break galaxies in the 
same field extended our knowledge of the X-ray properties of 
non-AGN dominated galaxies to z w 3 (Brandt et al. 2001a; 
Nandra et al. 2002). Lyman break systems, believed to be 
a combination of starbursts and low luminosity AGNs, arc 
found to have mean X-ray luminosities and Lx/Lb ratios 
comparable to those of the most luminous starbursts in the 
local universe (Brandt et al. 2001a). 



Although the studies above provide valuable informa- 
tion on the nature and the evolution of the X-ray emission 
in 'normal' galaxies, they focus on either local (m 20 Mpc) 
or higher redshift (z > 0.4) samples. Study of the X-ray 
properties of 'normal' galaxies in the intermediate redshift 
regime is still missing although critical in interpreting their 
evolutionary properties by providing a link between local 
and distant samples. In this paper we address this issue us- 
ing a shallow (2-10 ksec) XMM-Newton X-ray survey cover- 
ing an area of ~ 1.5 deg . These observations overlap with 
the 2dF Galaxy Redshift Survey (2dFGRS; mean redshift 
z ?s 0.1; Colless et al. 2001) providing a unique opportunity 
to study the mean X-ray properties of 'normal' galaxies at 
z ~ 0.1. Compared to previous studies, the 2dFGRS has the 
advantage of high quality spectra, accurate redshift mea- 
surements and reliable classifications into different galaxy 
types (E/S0, spirals) based on the observed spectroscopic 
properties. Moreover, sources with optical spectra revealing 
powerful AGNs can easily be discarded from the 'normal' 
galaxy sample. Our analysis of the mean X-ray properties 
of a well defined sample of 'normal' galaxies provides a ref- 
erence point at z « 0.1 for comparison with similar studies 
at higher redshifts. 

Section 2 presents the XMM/2dF survey, section 3 de- 
scribes the reduction of the X-ray data while section 4 details 
the selection criteria of the sample used in the present study. 
The stacking technique is outlined in section 5, the results 
are presented in section 6 and are discussed in section 7. 
Finally section 8 summarises our conclusions. Throughout 
this paper we adopt H Q = 65kms _1 Mpc -1 and q = 0.5. 



2 THE XMM/2DF SURVEY 

The North Galactic Pole F864 region (RA(J2000) = 13 h 41 m ; 
Dec.(J2000)=00°00') and the South Galactic Pole (SGP; 
RA(J2000)=00 h 57 m , Dec.(J2000) = -28°00') were surveyed 
by XMM-Newton between May and July 2002 as part of the 
Guaranteed Time program. The observations consist of 9 
pointings in the SGP and 7 pointings in F864 region each 
with an exposure time of w 2 — 10 ksec. The RA and DEC 
of the center of each pointing are listed in Table 1. The 
EPIC (European Photon Imaging Camera; Striider et al. 
2001; Turner et al. 2001) cameras were operated in full frame 
mode with the thin filter applied. 

Both the F864 and SGP regions overlap with the 
2dFGRst (Colless et al. 2001). The 2dFGRS is a large- 
scale spectroscopic survey that fully exploits the capabilities 
of the 2dF multi-fiber spectrograph on the the 4 m Anglo- 
Australian Telescope (A AT). This is an on-going project 
that aims to obtain high quality spectra and redshifts for 
250 000 galaxies brighter than an extinction corrected limit 
of bj < 19.4 mag over an area w 2000 deg 2 . Recently, data 
for about 100 000 galaxies have been released (100k data 
release). Part of this unprecedented spectroscopic database 
overlap with our XMM-Newton survey and hence optical 
spectroscopic information is available for galaxies in the 
F864 and SGP regions to the above magnitude limit. More- 
over, the central regions of both the F864 and SGP XMM- 

t http://msowww.anu.edu.au/2dFGRS/ 
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Newton survey have been observed by the ROSAT satellite 
(Griffiths et al. 1996; Georgantopoulos et al. 1996) 



3 XMM-NEWTON OBSERVATIONS 

The XXM-Newton data have been analysed using the Sci- 
ence Analysis Software (SAS 5.3). Event files for the PN 
and the two MOS detectors have been produced using the 
epchain and emchain tasks of SAS respectively. The event 
files were screened for high particle background periods by 
rejecting times with 0.5-10 keV count rates higher than 25 
and 15cts/s for the PN and the two MOS cameras respec- 
tively. These criteria excluded 5 fields (F864-4, SGP-1, 5, 

7 and 9) from the analysis that suffered by significantly el- 
evated and flaring particle background. The PN and MOS 
good time intervals for the remaining pointings are shown in 
Table 1. The differences between the PN and MOS exposure 
times are due to varying start and end times of individual 
observations. Only events corresponding to patterns 0-4 for 
the PN and 0-12 for two MOS cameras have been kept. To 
increase the signal-to-noise ratio and to reach fainter fluxes 
the PN and the MOS event files have been combined into a 
single event list using the merge task of SAS. 

Images in celestial coordinates with pixel size of 
4.35 arcsec have been extracted in the spectral bands 0.5- 
2keV (soft) and 2-8 keV (hard) for the merged event file. 
Exposure maps accounting for vignetting, CCD gaps and 
bad pixels have been constructed for each spectral band. 

Source detection is performed on the 0.5-8 keV image 
using the wavdetect task of CIAO with a false probability 
threshold of 10~ J (Freeman et al. 2002) corresponding to 
about 5a. A total of 508 X-ray sources have been detected 
to the limit /x(0.5 - 8keV) « lO^^ergs^ 1 cm" 2 . The 0.5- 

8 keV detection flux limit of individual XMM- Newton point- 
ings for a point source on-axis are shown in Table 1. We 
note that the false probability threshold of 10~ 5 for the to- 
tal band (0.5-8 keV) corresponds to a detection threshold 
of about 4cr in the soft band (0.5-2 keV). This limit trans- 
lates to a typical on-axis source flux of /x(0.5 — 2keV) « 

4 x 10~ 15 ergs -1 cm~ 2 in the soft 0.5-2 keV band. 

A byproduct of the source extraction algorithm is the 
construction of background maps for each spectral band. 
The mean background counts per pixel (estimated from the 
background maps) are also listed in Table 1. A detailed anal- 
ysis of the nature of the X-ray sources detected in the F864 
and SGP areas will be presented in a forthcoming paper 
(Georgantopoulos et al. in preparation). To convert counts 
to flux the Energy Conversion Factors (ECF) of individ- 
ual detectors are calculated assuming a power law spectrum 
with T = 1.7 and Galactic absorption N H = 2 x 10 20 cm~ 2 
appropriate for both the SGP and F864 fields. The mean 
ECF for the mosaic of all three detectors is estimated by 
weighting the ECFs of individual detectors by the respective 
exposure time. For the encircled energy correction, account- 
ing for the energy fraction outside the aperture within which 
source counts are accumulated, we adopt the calibration per- 
formed by Ghizzardi (2001a, 2001b). These studies use both 
PN and MOS observations of point sources to formulate the 
XMM-Newton PSF for different energies and off-axis angles. 
In particular, a King profile is fit to the data with parame- 
ters that are a function of both energy and off-axis angle. 



4 GALAXY SAMPLE 

The 'normal' galaxy sample used in the present study is com- 
piled from the 100k data release of the 2dF Galaxy Redshift 
Survey. A full description of the survey design, spectroscopic 
observations, data reduction and redshift measurements of 
the 2dFGRS dataset can be found at Colless et al. (2001). 
In brief, redshift estimation is performed using two indepen- 
dent methods: the first cross-correlates the spectra with a 
library of absorption-line galaxy and stellar templates after 
clipping emission lines, while the second finds and then fits 
emission lines. The automated redshift estimates are also 
confirmed by visual inspection. Occasionally, a manual red- 
shift is determined by fitting spectral features missed by 
the automated method (Colless et al. 2001). A quality flag, 
Q, measuring the reliability of the estimated redshifts has 
been assigned. Reliable redshifts have Q > 3, while Q — 2 
corresponds to a probable redshift and Q = 1 indicates no 
redshift measurement. 

Spectral classification of the 2dFGRS galaxies has 
been performed by applying Principal Component Analy- 
sis (PCA) techniques on the observed spectra (Folkes et al. 
1999; Madgwick et al. 2002). Spectral types are determined 
on the basis of the n parameter: 



r] — 0.5 x pa — pc2, (1) 

where pc\ , pc2 are the projections of the first two eigenspec- 
tra derived from PCA. Comparison of the n parameter of the 
2dFGRS galaxies with that derived for the observed spectra 
of the Kennicutt galaxy atlas (Kennicutt 1992) shows that 
there is a continuous relation between r) and morphological 
types (Madgwick et al. 2002, their Figure 4). We adopt the 
classification scheme of Madgwick et al. (2002) and group 
galaxies into four spectral types: E/S0, Sa, Sb and Scd. The 
adopted range of the r\ parameter for the above spectral 
types are presented in Table 2. "Unclassified" galaxies with 
no spectral classification are systems with low signal-to- 
noise spectra that did not allow spectral type determination. 

The present sample comprises a total of 235 2dFGRS 
galaxies with Q > 2 that lie within 13 arcmin off the center of 
the EPIC-PN camera to avoid the low signal-to-noise areas 
at the edge of the XMM-Newton field of view. Visual inspec- 
tion of the mosaiced images removes three galaxies close to 
CCD gaps. Moreover, two 2dFGRS sources (TGN335Z095 
and TGN 336Z241) lie close to an extended X-ray source, 
probably a cluster. These galaxies are not individually de- 
tected at X-ray wavelengths and are excluded from the 
stacking analysis described in section 5 to avoid contamina- 
tion of the stacked signal by emission not directly associated 
to the galaxies. A further 4 sources with stellar optical spec- 
tra are also excluded from the analysis reducing the final 
sample of 2dFGRS galaxies to 226. The distribution of the 
present sample of 2dFGRS galaxies to spectral types as well 
as the corresponding mean redshifts are shown in Table 2. 
We also note that the present sample is « 70 per cent com- 
plete. The incompleteness is due to Q < 2 galaxies as well 
instrumental limitations (e.g. fiber collisions). The optical 
B-band absolute magnitude distributions of ellipticala and 
spirals in the sample are shown in Figure 1. 
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Field Name 


RA 
(J2000) 


Dec 
(J2000) 


PN cxp. time 

(sec) 


MOS cxp. time 
(sec) 


/(0.5 - 8) a 


background 6 
(cnts) 


F864-1 


13 h 41 m 24.0 s 


+00°24'00" 

1 UU ^J'-t UU 


5779 


9974 


1.0 


0.33 


F864-4 


1 ^4111104 ns 


+oo°oo'oo" 

1 U U uu uu 










F864-5 


1 ^ h 4^ m nn n s 


+oo°oo'oo" 

1 u u uu uu 


1693 


4447 


2.8 


0.40 


F864-6 


1 ■jh^/moc ns 

-LiJ -L 1- UU.U 


+00°00'00" 

1 uu uu uu 


2766 


6493 


1.6 


0.29 


F864-7 


1 410104 ns 


— 00° 24' 00" 

UU iirt UU 


3459 


7139 


1.4 


0.27 


F864 8 


1 ^hiomn 4 ns 


— 00° 24' 00" 

— uu uu 


2109 


7276 


2.2 


0.62 


F864-9 


1 it/ii m i n n s 

1J Ttrt iJU.U 


— 00° 24' 00" 

UU iirt UU 




7725 


3.6 


0.46 


SGP-1 


uu <j<j ±ty.\j 


— 27°36'00" 

1 ou uu 










SGP-2 


00 h 57 m 00.0 s 


-27°36'00" 


4267 


8061 


1.1 


0.36 


SGP-3 


00 h 58 m 51.0 s 


-27°36'00" 


3181 


7729 


1.1 


0.20 


SGP-4 


00 h 55 m 19.0 s 


-28°00'00" 


4428 


8286 


1.1 


0.26 


SGP-5 


00 h 57 m 00.0 s 


-28°00'00" 










SGP-6 


00 h 58 m 51.0 s 


-28°00'00" 


6390 


10358 


1.0 


0.35 


SGP-7 


00 h 55 m 19.0 s 


-28°24'00" 










SGP-8 


00 h 57 m 00.0 s 


-28°24'00" 


7787 


3456 


0.9 


0.32 


SGP-9 


00 h 58 m 51.0 s 


-28°24'00" 











"Limiting 0.5-8 keV flux for detection (5ct) in units of 10 14 ergs 1 cm 2 

6 Mean 0.5-8 keV background counts per pixel from the background maps produced by the WAVDETECT task. 



Table 1. Observing Log of XMM-We«iion F864 and SGP survey 
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Figure 1. B-band absolute magnitude distribution of the spiral 
(dashed line) and the elliptical (continuous line) 2dFGRS sub- 
samples. For clarity the spiral galaxy histogram is shifted to the 
left by 0.1 magnitudes. 



Spectral 


V 


Number 


< z > 




type 




of galaxies 






all 




226 


0.110 


0.003 


E/S0 


r; < -1.4 


88 


0.105 


0.004 


Sa 


-1.4 < 77 < 1.1 


59 


0.104 


0.005 


Sb 


1.1 < ri < 3.5 


28 


0.092 


0.005 


Scd 


r; > 3.5 


19 


0.085 


0.007 


Unclassified 




32 


0.139 


0.013 



Table 2. Spectral classification scheme for the 2dFGRS galaxies 



5 STACKING PROCEDURE 

Stacking methods have been extensively used in X-ray as- 
tronomy to study the mean properties (e.g. flux, luminosity, 
hardness ratios) of well defined samples of sources that are 
otherwise too faint at X-rays to be individually detected 
(e.g. Nandra et al. 2002). 

In practice the X-ray counts (source-plus-background) 
at the position of each galaxy in the sample are added ex- 
cluding X-ray detected galaxies. The expected background 
contribution is estimated by summing the counts from re- 
gions around each source. Assuming Poisson statistics for 
the counts a significance level is estimated for the summed 
signal. 

To determine the size of the region within which the 
source-plus-background counts are extracted we adopt the 
empirical method described by Nandra et al. (2002). Tests 
are performed in which the radius of the circular aper- 
ture within which the source-plus-background counts are 
summed varies (from one trial to the next) between 8- 
18 arcsec. We find that a radius of 12 arcsec maximises the 
significance of the stacked signal and hence this optimal ex- 
traction radius is adopted for the analysis that follows. This 
extraction radius is about 3 times the on-axis HWHM of 
the XMM-Newton PSF (Hasinger et al. 2001). We also note 
that the visible disk of a Milky Way type galaxy at z=0.1 
subtends an angular radius of ~ 8 arcsec. 

To assess the significance of the stacked signal we esti- 
mate the background by averaging the counts in 60 12 arcsec 
apertures randomly positioned within annular regions cen- 
tered on the positions of the sample galaxies with inner and 
outer radii of 20 and 100 arcsec respectively. To avoid con- 
tamination from X-ray detections in the background esti- 
mation, regions close to X-ray sources (8 — 40 arcsec) are 
excluded. The smooth background map produced by the 
WAVDETECT task of CIAO was also used to estimate the 
background. Although no significant differences are found, 
the background estimated in the latter case is systematically 
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lower and hence, the detection significance of the stacked sig- 
nal is somewhat higher. We adopt a conservative approach 
and estimate the background from the science image rather 
than the smooth background map. However, it should be 
noted that use of the background map does not alter any of 
the conclusions. 

In practice, 2dFGRS galaxies associated (< 8arcsec) 
with X-ray sources detected in the total (0.5-8 keV) band 
above the false probability threshold of 1CF 5 are excluded 
from the stacking analysis. As already discussed in the pre- 
vious section this detection limit corresponds to a soft band 
(0.5-2 keV) detection threshold of about 4a which translates 
to a typical on-axis point source flux of /(0.5 — 2keV) ~ 
4 x 10~ 5 ergs - cm" . Such a low detection threshold is es- 
sential to study the mean properties of X-ray weak sources, 
the signal of which would otherwise be diluted by brighter 
ones. 

X-ray detected 2dFGRS galaxies may differ in their na- 
ture from X-ray faint galaxies. It is clear that summing 
counts of different classes of objects is meaningless. We find 
a total of 7 2dFGRS galaxies associated (< 8arcscc) with 
X-ray detections above the total band (0.5-8 keV) false prob- 
ability threshold of 10 . Two of them exhibit clear evidence 
for AGN activity such as broad emission line optical spec- 
tra, high X-ray-to-optical luminosity ratios (Lx/Lb > — 1) 
and X-ray luminosities Lx ~ 10 42 ergs _1 . Another three 
2dFGRS galaxies have narrow emission line optical spec- 
tra, Lx/Lb ~ — 1, relatively hard X-ray spectral prop- 
erties and soft band (0.5-2 keV) X-ray luminosities in the 
range 0.3 — 5 x 10 42 ergs _1 . The evidence above suggests 
obscured or weak AGN activity. Finally two 2dFGRS galax- 
ies have optical spectra exhibiting both absorption and nar- 
row emission lines, Lx/Lb ~ —2, soft X-ray spectra and 
L x ~ 2 x 10 41 ergs" 1 , suggesting either 'normal' galaxies 
or LLAGNs. Detailed analysis of these galaxies will be pre- 
sented in a forthcoming paper (Georgakakis ct al. in prepa- 
ration). 

Finally, counts associated with the wings of the PSF of 
bright X-ray sources might erroneously increase the stack- 
ing signal significance. Therefore, to avoid contamination 
of the stacking signal from nearby X-ray sources we have 
excluded from the analysis 2dFGRS galaxies that lie close 
(< 40arcsec) to X-ray detections (total of 21). There- 
fore, a total of 28 2dFRGS sources either lie close to (21; 
8 — 40 arcsec) or are associated with (total of 7; < 8 arcsec) 
X-ray detections. This reduces the number of galaxies used 
in the stacking analysis to 198. 



6 RESULTS 

The stacking results for the 2dFGRS galaxies are sum- 
marised in Table 3. The X-ray flux in this table is estimated 
from the raw net counts after correcting individual galax- 
ies for (i) the effect of vignetting estimated from the corre- 
sponding exposure map and (ii) the energy fraction outside 
the adopted aperture of 12 arcsec using the encircled energy 
corrections recently derived by Ghizzardi (2001a, 2001b) for 
both the PN and the MOS detectors (see section 3). The 
counts-to-flux conversion factor has been derived assuming 
a power-law with spectral index T — 1.7 and Galactic ab- 
sorption (see section 3). The X-ray luminosities in Table 3 



are estimated using the X-ray flux derived from the stacking 
analysis and the mean redshift of each subsample also listed 
in this table. 

The soft-band yields a statistically significant signal 
for both the whole sample and for the E/S0 and spiral 
galaxy sub-samples. On the contrary stacking of the hard 
band counts did not result in a significant signal for either 
of the samples. For the 2-8 keV band 3<r upper limits are 
estimated assuming Poisson statistics for the background 
stacked counts. Unfortunately, the estimated upper limits 
are not tight enough to strongly constrain the X-ray spec- 
tral properties of the spiral and elliptical galaxy subsamples. 

To further assess the confidence level of the stacked sig- 
nal in both the soft and the hard bands we perform extensive 
simulations: mock catalogs are constructed by randomising 
the positions of the 2dFGRS samples in Table 3 avoiding 
areas close to X-ray sources (< 40 arcsec). Each of the mock 
catalogs has the same number of sources as the 2dFGRS 
samples in Table 3. We then apply the stacking analysis 
technique outlined in section 5 to the mock catalogs and es- 
timate the significance of the stacked signal in the same way 
as in the real catalog. This is then repeated 10 000 times to 
get the distribution of the detection significances (as defined 
in section 5) for the mock catalogs. This provides an esti- 
mate of the probability of getting a statistically significant 
stacked signal by chance. The detection confidence levels 
estimated from the simulations are in excellent agreement 
with those shown in Table 3 based on Poisson statistics 

We also investigate the possibility that the observed 
stacked signal is primarily due to a few bright X-ray sources 
below the detection threshold. This is important especially 
if some of these sources are AGNs and hence, not representa- 
tive of the 'normal' galaxy population studied in this paper. 
Firstly, we find no evidence of broad emission lines in the 
optical spectra of the galaxies used in the stacking analysis. 
Therefore, the obtained stacked signal does not include any 
contribution from powerful AGNs/QSOs, although we can- 
not exclude the possibility of weak and/or heavily obscured 
AGNs present in the sample. We further explore this issue by 
assigning detection significances to individual 2dFGRS and 
then excluding sources above a given threshold. The detec- 
tion significance for each source is defined as the ratio of 
net source counts to the square root of the background. We 
note that this definition of detection significance is different 
to that used by the wavdetect task. 

For the whole sample we need to remove galaxies above 
2.7cr (total of 11) to obtain a low confidence level signal 
(2.6a). The 11 sources that appear to dominate the stack- 
ing signal have mean fx/ /b ~ —2 and Lx{0.5 — 2keV) « 
10 ergs - suggesting that the contribution from AGNs is 
small. Moreover, the optical spectra of the systems that ap- 
pear to dominate the stacking signal do not show any evi- 
dence for AGN activity (e.g. broad or high excitation lines) . 
The corresponding cutoffs for a low confidence stacked signal 
in the spiral and elliptical subsamples is 3<r. Deeper X-ray 
observations are required to further address this issue. 

Finally, we note that for the stacked signal we estimate 
an X-ray-to-optical flux ratio log fx/ fB ~ —2.5 typical of 
'normal' galaxies and Lx ~ 10 40 erg sec" 1 similar to that of 
Milky Way type galaxies. This also suggests that any contri- 
bution from AGNs in the stacking analysis is small. However, 
we note that the present data cannot be used to identify and 
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14.5 ±4.7 


< 21.4 


3.1 ± 1.0 


< 4.5 


8.5 
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87 
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3.3 


1.6 


2.5 ± 2.2 


< 27.8 


0.5 ±0.4 


< 2.2 


5.3 



a Source±Background counts in the soft (SB) and hard (HB) bands 
^Background counts in the soft (SB) and hard (HB) bands 

Significance of detection in background standard deviations in the soft (SB) and hard the (HB) bands 
d X-ray flux in units of 10~ 16 ergs -1 cm~ 2 in the soft (SB) and hard the (HB) bands 
e X-ray luminosity in units of 10 40 ergs _1 in the soft (SB) and hard the (HB) bands 
^B-band luminosity in units of 10 42 ergs -1 

Table 3. Stacking analysis results in the soft (0.5-2 keV) and hard (2-8 kcV) bands for different samples of 2dFGRS galaxies 
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A Shapley et al. 2002 

■ Hornschemeier et al. 2002 

• 2dFGRS Spirals 



0.4 0.6 0.B 

redshift 



Figure 2. log Lx/ Lb against redshift of spiral galaxies. The 
stacking analysis results for the 2dFGRS spirals are shown with 
the filled circle. The triangle is the median log Lx / Lg of local spi- 
rals (Shapley et al. 2001) selected to span an Lb range similar to 
that of the present sample. For this point the errorbars correspond 
to the 25th and 75th percentile of the log Lx/L B distribution of 
the Shapley et al. (2001) local spiral galaxy sample. The squares 
are the stacking analysis results for distant late type galaxies in 
the Chandra Deep Field North (squares; Hornschemeier 2002a). 



exclude from the stacking analysis weak or heavily obscured 
AGNs that do not manifest their presence in the optical 
galaxy spectra. For example, Comastri et al. (2002) have 
identified a class of X-ray bright optically inactive galaxies 
with optical spectra typical of 'normal' early type systems 
and X-ray properties indicative of obscured AGN activity. 
Our sample is likely to be contaminated by such systems. 
We note that because of their early type optical spectra 
sources like those found by Comastri et al. (2002) are likely 
to primarily affect the elliptical rather than the spiral galaxy 
subsample. 



7 DISCUSSION 

7.1 Optical— to— X-ray flux ratio 

Shapley, Fabbiano & Eskridge (2001) studied the X-ray 
properties of a large sample of local spirals using Einstein 



archival data. They performed multivariate statistical anal- 
ysis and report a non-linear correlation between optical and 
X-ray luminosities Lx ~ L B independent of spectral type. 
This contradicts previous studies, based on smaller samples, 
that found an almost linear Lx — Lb relation, interpreted 
as evidence for a link between X-rays and the evolved stellar 
population. To explain the non-linearity in Lx — L B Fab- 
biano & Shapley (2002) suggest two different processes af- 
fecting the X-ray emission of early and late type spirals: the 
presence of hot ISM similar to that of E/S0 in early types 
and dust obscuration in star-forming regions in late types. 

Brandt et al. (2001b) performed stacking analysis on 
a sample of 17 z £s 0.5 'normal' galaxies (both early and 
late type) using a 500 ks Chandra exposure of the HDF-N 
region. They estimate a mean 0.5-2 keV luminosity for the 
spiral galaxies in their sample of Lx = 8.5 x 10 39 erg s _1 con- 
sistent with the calculations in the present paper. This Lx 
is elevated compared to local spirals but this is attributed to 
the increased mean Lb in their sample rather than evolution. 
Similarly, Hornschemeier et al. (2002a) used stacking analy- 
sis techniques to study the X-ray properties of the spiral and 
irregular type systems in the Chandra Deep Field North at 
redshifts z 0.4 — 1.5. They argue that the mean Lx I Lb in 
their sample is elevated compared to that of Shapley et al. 
(2001) (selected to span the same Lb range) suggesting evo- 
lution. However, the uncertainties in their Lx/ Lb estimates 
do not allow firm conclusions to be drawn. 

The present sample with well defined selection crite- 
ria is intermediate in redshift (z ~ 0.1) to that of Shapley 
et al. (2001; z « 0.004) and Hornschemeier et al. (2002a; 
z > 0.4), providing a low- z reference point for comparison 
with higher-z studies. Figure 2 plots the log Lx / Lb against 
redshift comparing the results from the studies above. For 
the Shapley et al. (2001) spirals the median Lx / Lb is es- 
timated using galaxies over the same range of Lb as the 
present spiral galaxy subsample. Within the uncertainties 
there is fair agreement in the Lx/Lb of local and z m 0.1 
2dFGRS spirals. 

Moreover, it is clear from Figure 2 that the Lx/Lb 
of z ~ 0.1 spirals in the present study are also consistent, 
within the errors, with those of Hornschemeier et al. (2002a) 
at z = 0.4 — 1. Although the mean Lx of their spirals in- 
creases with redshift from z ~ 0.4 to z ~ 1 this is com- 
pensated by a similar increase in the mean Lb- At z > 1 
Hornschemeier et al. (2002a) find that both Lx/Lb and 
Lx are elevated compared to our estimates by a factor of 
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4 suggesting evolution. Although the large errors hamper a 
secure interpretation, it should be noted that this increase, 
if real, cannot be attributed to differences in Lb between our 
sample and that of Hornschemeier et al. (2002a) at z > 1. 
Therefore, our analysis combined with the results of Horn- 
schemeier et al. (2002a) suggest little or no evolution of the 
spiral galaxy Lx relative to Lb at least to z « 1. 

7.2 X-ray emissivity and star- formation rate 

The emissivity of galaxy samples selected at various wave- 
lengths (UV, optical, radio) over a range of redshifts is one 
of the most powerful and widely used evolutionary tests (e.g. 
Lilly et al. 1996; Connolly et al. 1997). The X-ray emissiv- 
ity, jx , of different galaxy types (including non-AGN domi- 
nated systems) has been estimated locally by Georgantopou- 
los, Basilakos & Plionis (1999). However, the evolution of 
jx at higher redshifts is still poorly constrained mainly be- 
cause 'normal' galaxies are faint at X-ray wavelengths and 
their detection, especially at high redshift, is difficult even 
with the improved sensitivities of the XMM-Newton and the 
Chandra observatories. 

The 2dFGRS sample at z ft* 0.1 can provide, for the first 
time, a non-local estimate of the X-ray emissivity of 'normal' 
galaxies. For the jx determination, we first need to estimate 
the effective cosmological volume probed by the 2dFGRS 
survey corrected for optical magnitude incompleteness (e.g. 
bj < 19.4 mag) using the standard 1/V ma x formalism (e.g. 
Lilly et al. 1996). Corrections have also been applied tak- 
ing into account incompleteness of the 2dFGRS sample due 
to both poor S/N spectra (Q < 3, unclassified objects) at 
faint magnitudes and instrumental limitations (e.g. fiber col- 
lisions) that did not allow a fraction of the target galaxies to 
be observed. We further assume that all 2dFGRS galaxies 
of a given spectral type have the same mean Lx/Lb ratio 
estimated in section 6. 

For 2dFGRS spirals we estimate a soft-band 0.5-2 keV 
emissivity of (0.4 ± 0.3) x 10 38 in fair agreement with that 
of local HII galaxies, (0.50 ± 0.06) x lO^ergsec" 1 Mpc~ 3 , 
estimated by Georgantopoulos et al. (1999; scaled to the 
0.5-2 keV band). For the elliptical galaxy sub-sample we find 
jx = (1.2 ± 0.4) x lO^ergsec^ 1 Mpc~ 3 . Although the un- 
certainties are large this is higher than the jx of local pas- 
sive galaxies, (0.53± 0.06) x 10 38 erg sec" 1 Mpc~ 3 , estimated 
by Georgantopoulos et al. (1999). These authors convolved 
the local optical luminosity function of the Ho, Filippenko 
& Sargent (1995) sample with the corresponding Lx — Lb 
relation based on Einstein data. The Ho et al. (1995) sam- 
ple has the advantage of high S/N nuclear optical spectra 
and hence, reliable spectral classifications (Ho, Filippenko & 
Sargent 1997) necessary to study the properties of non-AGN 
galaxies. 

Assuming the same mean 0.5-2 keV X-ray flux (esti- 
mated from stacking; see Table 1) for all 2dFGRS spi- 
ral galaxies and an incompleteness of ps 70 per cent 
(see section 3) we estimate their XRB contribution to be 
ps 0.4 per cent. This calculation assumes that the spec- 
trum of the diffuse X-ray background follows the form 
9 x S"°' 4 keVcm" 2 s~ 1 sr" 1 keV" 1 (Gendreau et al. 1995). 
This small XRB fraction produced by 'normal' spirals at 
z ps 0.1 is broadly consistent with studies of 'normal' galax- 
ies at higher redshifts using deep X-ray surveys. Horn- 



schemeier et al. (2002a) argue that they have identified ps 1 
per cent of the soft XRB as arising from spirals not indi- 
vidually detected in deep Chandra surveys. Alexander et 
al. (2002) studied X-ray detected starbursts in the redshift 
range 0.4-1.3 and argued that they produce ps 2 per cent of 
the 0.5-8 keV XRB. 

In addition to the above observational constraints, a 
number of studies attempt to predict the total XRB frac- 
tion produced by an evolving population of star-forming 
galaxies. Ranalli et al. (2003) used a nearby well defined 
sample of star-forming galaxies with available ASCA and 
BeppoSAX data to derive an almost linear relation between 
X-ray luminosity and radio power. Using this correlation 
and the radio source counts at sub-mjy flux densities (be- 
lieved to be dominated by starbursts) they predict the con- 
tribution of these systems to the XRB. Integrating their 
predicted soft band (0.5-2 keV) X-ray source counts we es- 
timate a 0.5-2 keV XRB fraction produced by starbursts in 
the range ps 10 — 15 per cent depending on the faint flux 
limit of the integration. Persic & Rephaeli (2003) carefully 
modeled the X-ray spectral energy distribution of nearby 
star- forming galaxies to estimate a 2-10 keV XRB contribu- 
tion of 5 — 11 per cent depending on the assumed evolution. 
Although they do not estimate the soft band XRB contribu- 
tion the fractions above will increase at softer energies. Grif- 
fiths & Padovani (1990) converted the 60 — /im luminosity 
function of IRAS star-forming galaxies to X-ray wavelengths 
and argued that star-forming galaxies with mild evolution 
could produce 20-30 per cent of the XRB. Similarly, Pearson 
et al. (1997) used IRAS star-forming galaxies to fit the X-ray 
counts at faint flux limits and estimated that such systems 
could contribute as much as 30-50 per cent of the XRB. The 
studies above clearly demonstrate that there is significant 
uncertainty in the XRB fraction produced by star-forming 
galaxies critically depending on the adopted evolution. 

The present study of spiral galaxies to z ps 0.1 can- 
not on its own be used to either constrain the evolution of 
these systems with z or estimate their total XRB contri- 
bution extrapolated to high z. However, comparison of our 
results at z ps 0.1 with those obtained at higher redshifts 
(e.g. Hornschemeier et al. 2002) indicate X-ray evolution for 
these systems. Indeed, the almost constant Lx/Lb versus 
z relation in Figure 2 can be interpreted as evidence that 
Lx evolves in at least the same rate as Lb - It is an interest- 
ing exercise to attempt to constrain the X-ray evolution of 
star-forming galaxies using the X-ray emissivity of 2dFGRS 
spirals at z ps 0.1 and assuming that the bulk of the soft 
XRB that remains unaccounted for in deep X-ray surveys is 
produced by star-forming processes. Deep ROSAT surveys 
(Hasinger at al. 1998) to the limit ps 10 _15 ergs _1 cm~ 2 have 
resolved 68-81 per cent of the soft (< 2 keV) XRB into dis- 
crete sources. Most of them are AGNs (78 per cent; Schmidt 
et al. 1998) with a small fraction being galaxy groups or clus- 
ters (6 per cent; Lehmann et al. 2001). More recently, deep 
Chandra surveys (e.g. Mushotzky et al. 2000; Brandt et al. 
2001c; Giacconi et al. 2002) extended the study of the X-ray 
source population to the limits ps 10~ 17 — 10 _16 crgs _1 cm~ 2 . 
For example, Mushotzky et al. (2000) resolved 6-13 per cent 
of the soft band XRB to sources also believed to be AGNs 
with fluxes in the range ps 10~ 16 — 10 _15 ergs _1 cm~ 2 . Only 
a small fraction of the XRB (1-2 per cent) is believed to 
be produced by star-forming galaxies at the flux limits of 
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the ultra-deep Chandra surveys (Hornschemeir et al. 2002a; 
Alexander et al. 2002). The evidence above suggests that 
74-96 per cent of the soft band XRB has been resolved and 
is produced by AGNs, groups or clusters allowing a maxi- 
mum of 6-26 per cent to star-forming galaxies. We employ 
the spiral galaxy jx at z « 0.1 to constrain the maximum 
allowed evolution for these systems. Under the assumption 
that the X-ray emissivity follows the form (l + z) k to a max- 
imum formation redshift z m ax = 4, we find that the index k 
should be < 1.5 and < 3.2 for the XRB fraction produced 
by these systems to be < 6 and < 26 per cent respectively. 
We note that surveys at optical wavelengths estimate k ^ 3 
(e.g. Lilly et al. 1996). 

We stress that the XRB fraction due to star-forming 
galaxies quoted above (and hence, the derived evolution pa- 
rameter k) is an upper limit. Although, a number of studies 
suggest that as much as 20-50 per cent of the XRB could 
be produced by starbursts (e.g. Griffiths & Padovani 1990; 
Pearson et al. 1997) more recent estimates suggest modest 
soft band fractions (e.g. Persic & Rephaeli 2003; Ranalli 
et al. 2003). Deep wide field surveys selecting star- forming 
galaxies at higher redshifts using homogeneous and well de- 
fined criteria are required to constrain the X-ray evolution of 
these systems and hence, their XRB contribution. This pa- 
per provides such a sample at low redshifts for comparison 
with future studies. 

Finally, the X-ray emission of spirals not dominated 
by AGN is associated to the star-formation activity of the 
parent galaxy (Fabbiano & Shapley 2002; Read & Ponman 
2001). Consequently, the X-ray emission of these systems 
can be used to probe both the individual and the global 
star- formation rate (SFR). However, the quantitative link 
between X-ray luminosity and SFR is rather tenuous, de- 
pending on poorly constrained parameters such as the Ini- 
tial Mass Function and the nature of the stellar sources 
responsible for the X-ray emission (e.g. Ghosh & White 
2001). More theoretical and observational work is clearly 
needed in this field. Here, we adopt the empirical 0.5-2 keV 
X-ray luminosity to SFR conversion derived by Ranalli et 
al. (2003). For the X-ray luminosity of the spirals in the 
present sample we find a mean SFR of ~ 1.3M0yr _1 
while the star-formation density of the present sample at 
z « 0.1 is « 0.009 M yr _1 Mpc~ 3 . This is lower than pre- 
vious estimates using galaxy samples selected at different 
wavelengths. However, given the large uncertainties involved 
in converting X-ray luminosity to star-formation rate the 
agreement is surprising. This suggests that deep X-ray sur- 
veys detecting an increasing number of 'normal' galaxies can 
be used as probes of the SFR of individual systems. 



8 CONCLUSIONS 

In this paper we study the mean X-ray properties of 'nor- 
mal' galaxies using stacking analysis techniques. We use a 
shallow (2-10 ksec) XMM-Newton survey covering an area of 
w 1.5 deg 2 and overlapping with the 2dF Galaxy Redshift 
Survey. The 2dFGRS sample has the advantage of high qual- 
ity spectra and hence, accurate redshifts and reliable classi- 
fications into early-type and late-type galaxies. 

In the stacking analysis we use a sample of « 200 2dF- 
GRS galaxies with a mean redshift of z w 0.1. In the soft 



0.5-2 keV band we detect a statistically significant signal for 
both the whole sample (~ 6 <r) and the early (ss 4<r) and late 
(w 3 a) type galaxy sub-samples. On the contrary, no signal 
is detected for any of the above samples in the hard 2-8 keV 
band. The main conclusions from this study are summarised 
below: 

• We find that the mean Lx /Lb ratio of the present sam- 
ple at z « 0.1 is consistent within the errors with both local 
(< 100 Mpc) and distant (z w 1) samples after accounting 
for differences in the mean Lb of the samples. This suggests 
little or no evolution of the X-ray emission mechanisms rel- 
ative to the optical at least to z « 1. 

• The mean X-ray emissivity, jx , of 'normal' galaxies at 
z « 0.1 is found to be consistent with that of local HII 
samples. We also find that the present sample of 2dFGRS 
spirals contributes about 0.4 per cent to the soft band XRB. 
This low fraction is in broad agreement with previous stud- 
ies estimating the XRB produced by star-forming spirals at 
moderate redshifts. 

• For star-forming galaxies assuming an X-ray luminos- 
ity evolution of the form (1 + z) k to a maximum redshift 
Zmax = 4 we find k < 3. The k values are constrained by the 
relatively low fraction of the soft X-ray background that re- 
mains unresolved by deep surveys (6-26%). Higher k values 
will result in an overproduction of the soft X-ray background 
given the fraction already attributed to AGN, groups and 
clusters. 

• Using the mean X-ray properties of the spiral galaxy 
sub-sample and assuming an empirical X-ray luminosity to 
SFR conversion factor we estimate a global SFR density of 
0.009±0.007 M yr _1 at z « 0.1. Although the uncertainties 
are large this is lower than previous results based on galaxy 
samples selected at different wavelengths. 
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